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ABSTRACT 

We study resonant leptogenesis at the electroweak phase transition in the singlet Majoron 
model with right-handed neutrinos. We consider a scenario, where the SM gauge group 
and the lepton number break down spontaneously during a second-order electroweak phase 
transition. We calculate the flavour- and temperature-dependent leptonic asymmetries, by 
including the novel contributions from the transverse polarisations of the and Z bosons. 
The required resummation of the gauge- dependent off-shell heavy-neutrino self-energies is 
consistently treated within the gauge-invariant framework of the Pinch Technique. Taking 
into consideration the freeze-out dynamics of sphalerons, we delineate the parameter space 
of the model that is compatible with successful electroweak resonant leptogenesis. The 
phenomenological and astrophysical implications of the model are discussed. 

PACS numbers: ll.30.Er, 14.60.St, 98.80.Cq 
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1 Introduction 



Leptogenesis [1] provides an elegant framework to consistently address the observed Baryon 
Asymmetry in the Universe (BAU) [2] in minimal extensions of the Standard Model (SM) [3]. 
According to the standard paradigm of leptogenesis, there exist heavy Majorana neutrinos 
of masses close to the Grand Unified Theory (GUT) scale Afcux ~ 10"*^^ that decay out 
of equilibrium and create a net excess of lepton number (L), which gets reprocessed into 
the observed baryon number (B), through the (i? + L)-violating sphaleron interactions [4]. 
The attractive feature of such a scenario is that the GUT-scale heavy Majorana neutrinos 
could also explain the observed smallness in mass of the SM light neutrinos by means of 
the so-called seesaw mechanism [5]. 

The original GUT-scale leptogenesis scenario, however, runs into certain difficulties, 
when one attempts to explain the flatness of the Universe and other cosmological data [2] 
within supergravity models of inflation. To avoid overproduction of gravitinos G whose late 
decays may ruin the successful predictions of Big Bang Nucleosynthesis (BBN), the reheat 
temperature Trch of the Universe should be lower than 10^-10*' GeV, for = 8-0.2 TeV [6]. 
This implies that the heavy Majorana neutrinos should accordingly have masses as low as 
Treh ~ 10^ GeV, thereby rendering the relation of these particles with GUT-scale physics 
less natural. On the other hand, it proves very difficult to directly probe the heavy-neutrino 
sector of such a model at high-energy colliders, e.g. at the LHC or ILC, or in any other 
foreseeable experiment. 

A potentially interesting solution to the above problems may be obtained within the 
framework of resonant leptogenesis (RL) [7]. The key aspect of RL is that self-energy ef- 
fects dominate the leptonic asymmetries [8], when two heavy Majorana neutrinos happen 
to have a small mass difference with respect to their actual masses. If this mass difference 
becomes comparable to the heavy neutrino widths, a resonant enhancement of the leptonic 
asymmetries takes place that may reach values 0{1) [7,9]. An indispensable feature of RL 
models is that flavour effects due to the light-to-heavy neutrino Yukawa couplings [10] play 
a dramatic role and can modify the predictions for the BAU by many orders of magni- 
tude [11,12]. Most importantly, these flavour effects enable the modelling [12] of minimal 
RL scenarios with electroweak-scale heavy Majorana neutrinos that could be tested at the 
LHC [13, 14] and in other non- accelerator experiments, while maintaining agreement with 
the low-energy neutrino data. Many variants of RL have been proposed in the litera- 
ture [15,16], including soft leptogenesis [17] and radiative leptogenesis [18]. 

In spite of the many existing studies, leptogenesis models face in general a serious 
restriction concerning the origin of the required CP and L violation. If CP or L violation 
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were due to the spontaneous symmetry breaking (SSB) of the SM gauge group, a net 
L asymmetry could only be generated during the electroweak phase transition (EWPT), 
provided the heavy Majorana neutrinos are not too heavy such that they have not already 
decayed away while the Universe was expanding]*] 

In this paper we show how RL constitutes an interesting alternative to provide a 
viable solution to the above problem as well. For definiteness, we consider a minimal 
extension of the SM with right-handed neutrinos and a complex singlet field S. The 
model possesses a global lepton symmetry U(l)i which gets spontaneously broken through 
the vacuum expectation value (VEV) of S, giving rise to the usual AL = 2 Majorana 
masses. Because of the SSB of the U(l);, the model predicts a true massless Goldstone 
boson, the Majoron. Therefore, this scenario is called the singlet Majoron model in the 
literature [20,21]. Depending on the particular structure of the Higgs potential, the VEV 
of S may be related to the VEV of the SM Higgs doublet $. Such a relation, for example, 
arises if the bilinear operator S*E is small or absent from the Higgs potential. In this case, 
the breaking of L occurs during the EWPT. For the model under study and given the LEP 
limit [22] on the SM Higgs boson Mj^ > 115 GeV, the EWPT is expected to be second 
order and hence continuous from the symmetric phase to the broken one [23] . 

We should now notice that all SM fermions and right-handed neutrinos have no 
chiral masses above the EWPT and therefore the generation of a net leptonic asymmetry 
is not possible. Consequently, in this model successful baryogenesis can result from RL at 
the EWPT. Although the singlet Majoron model that we will be studying here violates 
CP explicitly, the results of our analysis can straightforwardly apply to models with an 
extended Higgs sector that realise spontaneous CP violation at the electroweak scale. 

The paper is organised as follows: Section [2] presents the basic features of the singlet 
Majoron model with right-handed neutrinos, including the interaction Lagrangians that are 
relevant to the calculation of the leptonic asymmetries in Section [31 Moreover, in Section [3] 
we consider the novel contributions to the leptonic asymmetries, coming from the transverse 
polarisations of the and Z bosons. In the same context, the resummation of the 
gauge-dependent off-shell heavy-neutrino self-energies [24, 25] (which remains an essential 
operation in RL) is performed within the so-called Pinch Technique (PT) framework [26]. 
In Section m we analyse the Boltzmann dynamics of the sphaleron effects on RL and present 
predictions for the BAU. Section [5] is devoted to the phenomenological and astrophysical 
implications of the singlet Majoron model. Finally, Section [6] contains our conclusions. 

*An exception to this argument may result from a phase transition that is strongly first order. However, 
such a scenario is not feasible within the SM with singlet neutrinos [19] (see also our discussion below). 
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2 The Singlet Majoron Model 



Here we describe the basic features of the singlet Majoron model [20,21] augmented with 
a number n^ of right-handed neutrinos (with a = 1,2, .. . ,n/{) that will be relevant 
to our study. As mentioned in the introduction, the singlet Majoron model contains one 
complex singlet field S in addition to the SM Higgs doublet $. Although E is not charged 
under the SM gauge group 811(2)^ ® U(l)y, it still carries a non-zero quantum number 
under the global lepton symmetry U{l)i- More explicitly, the scalar potential of the model 
is given by 

- Cv = m|<l>^<l> + m|E*E + ^ ('^^*)' + Y " <5$^$S*S. (2.1) 

In order to minimise the potential fl2.ll) . we first linearly decompose the scalar fields as 
follows: 

/ G+ 



Then, the extremal or tadpole conditions may easily be calculated by 

If m| or m| are negative, the tadpole conditions (12.31) and (12.41) imply that the ground 
state of the scalar potential breaks spontaneously the local 811(2)/, (8)U(l)y and the global 
U{l)i symmetries, through the non-zero VEVs v and w, respectively. 

Expanding the fields $ and S about their VEVs, we obtain three would-be Goldstone 
bosons and G^, which become the longitudinal polarisations of and Z bosons, 
and one true massless Goldstone boson J associated with the SSB of U(l);. This massless 
CP-odd field J is called the Majoron in the literature [20,21]. In addition, there are two 
CP-even Higgs fields H and S, whose masses are determined by the diagonalisation of the 
mass matrix 

= ( Y ) , (2.5) 

where Ai^ is defined in the weak basis (0, a). The Higgs mass eigenstates H and S are 
related to the states (p and a, through the orthogonal transformation: 

. . (2.6) 

a 




with = Sp/cp = v/w and 

t20 = f-^ . (2.7) 

In the above we used the short-hand notation: = sin a;, c^. = cos a; and t^. = tanx. 
Moreover, the squared mass eigenvalues of the CP-even H and S bosons may easily be 
calculated from in (12.51) and are given by 

,2 



Mis = y 



The requirement that Mfj g be positive gives rise to the inequality conditions, 



(2.8) 



A$,s > , A$As > 5% (2.9) 

for the quartic couplings of the potential. In this context, we note that if |m|| <^ 

(5/A$) |m|| such that m| can be completely neglected in the scalar potential, the VEV w 

of S is then entirely determined by the VEV f of $ and the quartic couplings As and 6, 

viz. ^ 
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w = J— V . (2.10) 



This is an interesting scenario, since the ratio tp = v/w = y Xt,/6 does not strongly depend 
on the temperature T, as opposed to what happens to the VEVs v and w individually. In 
fact, as long as A$^s, ^ <^ 1, the thermally-corrected effective potential can be expanded, to 
a very good approximation, in powers of T^/m|. In such a high-T expansion, the quartic 
couplings of Cy turn out to be T-independent [27] and hence tp does not depend on T. 

We now turn our attention to the neutrino Yukawa sector of the model, which is 
non-standard. After SSB, it is given in the unitary gauge by 

- Cy = - l'iL{mD)ial^aR + ^ t ^'^ ^aR {'n^M)aP ^PR + H.C., (2.11) 

V 2w 
where summation over repeated indices is understood. Hereafter we use Latin indices 
to label the left-handed neutrinos, e.g. un, and Greek indices for the right-handed ones, 
e.g. i^aR- Observe that the spontaneous breaking of U(l)/ generates lepton-number- violating 
AL = 2 Majorana masses {mM)ap in addition to the lepton-number-preserving AL = 
Dirac masses {m£))ia. 

The model under discussion predicts a number (3 + n^j) of Majorana neutrinos which 
we collectively denote by nj, with I = i , a. Their physical masses are obtained from the 
diagonalisation of the neutrino mass matrix 

= \ °, 1 , (2.12) 
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by means of the unitary transformation U'^'^ M." U" = Ai", where tVI'^ is a non- negative 
diagonal matrix. The neutrino mass eigenstates (ri/)/? and (n/)^ are related to the states 
i^iL, (i^iL)'^, I'aR and {faR)^ through 

j = U'ij {nj)n , [ ] = ^iJ Ml • (2.13) 

Assuming the seesaw hierarchy {mD)ia/ {^M)ai3 1, the model predicts 3 light states 
that are identified with the observed light neutrinos {rii = z/j), and a number n^j of heavy 
Majorana neutrinos (n^ = A^^) with masses of order {rnM)ap = Pa/dU), where = Ppa are 
the Yukawa couplings of S to right-handed neutrinos. 

To obtain an accurate light and heavy neutrino mass spectrum within the context of 
models of electroweak RL, it is important to go beyond the leading seesaw approximation. 
To this end, we need first to perform a block diagonalisation and cast M.^ into the form: 

( T ) ■ '-^) 

This can be achieved by introducing the unitary matrix V [28]: 

[ -eih + ce)-'/' (In, + en-'/' ) ' ^' ' 

where ^ is an arbitrary 3 x n^j matrix. The expressions (13 + ^*^-^)"^/^ and (In, +^"^^*)^^''^ 
are defined in terms of a Taylor series expansion about the A/" x A/" identity matrix l_\f. 
These infinite series converge provided the norm ||^|| is much smaller than 1, where ||^|| = 



Tr(,^,^"l"). This condition is naturally fulfilled within the seesaw framework [29]. Block 
diagonalisation of the matrix Ai'^ given in f l2.12p implies that the {12} block element of 
V^Ai'^V vanishes, or equivalently that 

mo - iniM - imlC = . (2.16) 

Equation fl2.16p determines ^ in terms of m£, and m^f. It can be solved iteratively, with 
the first iteration given by 

^ = moml} - mornlj m^m*jjmlY' ml^ . (2.17) 

Note that the second term on the RHS of fl2.17p is suppressed by the ratio of the light- 
to-heavy neutrino masses and can thus be safely neglected in numerical estimates. Upon 
block diagonalisation, the block mass "eigen-matrices" are 

= (l„^ + e^y'^' {rriM + m^C + ^^run) (in, + eC)''^' , (2.18) 



m'^ = - (I3 + a^)" {rnoe + irn^ - iniMe) (U + T^' 

= -em^e^, (2.19) 
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where we used (I2.16P to arrive at the last equahty of fl2.19p . Keeping the leading order 
terms in an expansion of in powers of m/jm^/, we find that 



= niM + ^ (^m^D^M^D + m^m^^m*^ j , 



m 



-1 N -1 T 



(2.20) 

These last expressions are used to calculate the light and heavy neutrino mass spectra of 
the RL scenarios discussed in Section HI 

In order to calculate the leptonic asymmetries in the next section, we need to know the 
Lagrangians that govern the interactions of the Majorana neutrinos rii and charged leptons 
I = e, fi, T with: (i) the and Z bosons; (ii) their respective would-be Goldstone bosons 
and G; {Hi) the CP-odd Majoron particle J; (iv) the CP-even Higgs fields H and 5*. 
In detail, these interaction Lagrangians are given by [21] 



V2 



W-^lBin^PLUi + H.c, 



J^"" G' IBu (rriiPL - rrii Pr ) rij + H.c. 

V 2 My/ ^ ' 



4 cos Q, 
9 



G 



4M, 



G rij 



w 



Gij (mi Pl - mj Pr) + C}j (mj Pl - rm Pr 



(2.21) 
(2.22) 
(2.23) 

nj , (2.24) 



^9^ 



-I- 5ij m/75 



tg J rii 



Cij (mi Pl - mj Pr) + C}j {nij Pl - rm Pr 



nj 



(2.25) 



C 



9v 



H 



4M, 



w 



Cij (mi Pl + nij Pr) + C}j (m j Pl + mj Pr 



5lj mi 75 



(2.26) 



9v 



+ 



te + tf3 



se + cetp) S rii 
Sij mi 75 



Cij (mi Pl + mj Pr) + C*ij (mj Pl + mj Pr 



nj 



where Pl,r = | (1 =F 75), 9w is the SU(2)i gauge coupling of the SM and 



Bn = viur 



'Ik ^kl 5 



Cij - U^iUj^} 



(2.27) 



(2.28) 
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In (12.281) V'' is a 3-by-3 unitary matrix that occurs in the diagonahsation of the charged 
lepton mass matrix AiK Without loss of generahty, we assume throughout the present 
study that A^' is positive and diagonal, which implies that V'' = I3. Finally, we comment 
on the limit of — > 0. It is easy to see from (12. 7p that this limit leads to to and the 
fields S and J decouple from matter; only the Higgs field H couples to Majorana neutrinos 
and to the rest of the SM fermions (cf. [13]). 

3 Leptonic Asymmetries 

In this section we calculate the leptonic asymmetries produced by the decays of the heavy 
Majorana neutrinos during a second-order EWPT. The novel aspect of such a calculation 
is that, in stark contrast to the conventional leptogenesis scenario, the and Z bosons 
also contribute to the decays and leptonic asymmetries of the heavy Majorana neutrinos. 
This fact raises new issues related to the gauge invariance of off-shell Green functions which 
are here addressed within the so-called Pinch Technique (PT) framework [26]. 

Since sphalerons act on the left-handed SM fermions converting an excess in leptons 
into that of baryons, we only need to consider the decays of the heavy Majorana neutrinos 
Na into the left-handed charged leptons and light neutrinos uil. In detail, we have 
to calculate the partial decay width of the heavy Majorana neutrino into a particular 
lepton flavour /, 

r'^^ = r(iv, ^ /- w+, G+) + r(iv, ^ z, g, j, h, s) . (3.1) 

To compute F^y^, it proves more convenient to first calculate the absorptive part T,^^{^) 
of the heavy Majorana-neutrino self-energy transition Nj^ — > in the Feynman-'t Hooft 
gauge ^ = 1, where is the 4-momentum carried by Na^/s- The Feynman-'t Hooft gauge is 
not a simple choice of gauge, but the result obtained in the gauge-independent framework 
of the PT [26], within which issues of analyticity, unitarity and CPT invariance can self- 
consistently be addressed [24,25]. 

Neglecting the small charged-lepton and light-neutrino masses, acquires the 

simple spinorial structure: 

E5(^) = A^p{s) ^Pl + Kf^is) ^Pr , (3.2) 

where s = is the squared Lorentz-invariant mass associated to the self-energy transi- 
tion A^^ —>■ Na. Considering the Feynman graphs shown in Fig. [1] and the interaction 
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ll , VlL 

(a) 



Nr. 



(b) 



Figure 1: Feynman graphs that determine the 1-loop absorptive part ^^^(s) of the heavy 
Majorana-neutrino self-energy Sq,^(^). 



Lagrangians fl2.2ip - fl2.27l) . the absorptive transition amplitudes Aai3{s) are calculated to 
be 



+ 



32 



l=e,fj. 

rriN^ rriN^ 



4 1 ^ e s- Ml) + 1 



Ml 



(s - Mi 



Mir 



M^ \ M^ 



0(s - M^) + 1 - 



+ (C0 - setpf 1 



^\6{s-Ml) + (., + c,t,ffl-^V^(s-M|; 



(3.3) 



where = 5f^/(47r) is the SU(2)i fine-structure constant and 6{x) is the usual step 
function: 9{x) = 1 for a; > 0, whilst 6{x) = if x < 0. In the calculation of ^^^(s), 
we used the fact that Bia = C^ia + ^iCl^^), which is an excellent approximation in the 
physical charged-lepton mass basis. 

We should bear in mind that all masses involved on the RHS of (13. 3 p depend on 
the temperature T, through the T-dependent VEVs v{T) and w{T) related to the Higgs 
doublet $ and the complex singlet S, respectively [cf. (14.80 and fl4.17p ]. In the symmetric 
phase of the theory, i.e. for temperatures above the electroweak phase transition, these 
VEVs vanish and the absorptive transition amplitude becomes 



(m^ my)^p ( . , t} 



(3.4) 



Note that this last formula is only valid in the weak basis in which the Majorana mass 
matrix niM is diagonal. 

To account for unstable-particle-mixing effects between heavy Majorana neutrinos, 
we follow [7, 9] and define the resummed effective couplings Bia and their CP-conjugate 
ones B'I^ related to the vertices W'liNa and W^{Il)^ Na, respectively. For a symmetric 
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model with 3 left-handed and 3 right-handed neutrinos, the effective couplings Bia exhibit 
the same analytic dependence on the absorptive transition amplitudes Aa/B as the one found 
in [9]|l| 

3 

Bla = Bia — "i X! Bijj (3.5) 

/3,7=1 

X 







+ mpAp^{maA^a + m^Aa^) 


m?^ - TTLp + 2i mlApf3 + 2i ImRa^ (ml 


A/3^ 


2 + mpm-yKeA^p^^ 





where all transition amplitudes Aap, Ap^ etc are evaluated at s = rrij^^ = and 



2 

TTl 

Moreover, l^a/j-yl is the modulus of the usual Levi-Civita anti-symmetric tensor. The re- 
spective CP-conjugate effective couplings B^^ are easily obtained from (13.51) by replacing 
the ordinary iy~-boson couplings Bia and Aap{s) by their complex conjugates. In the 
decoupling limit of m^rg ^ m^^^^, we recover the analytic results known for a model with 
2 right-handed neutrinos [7,9], where the effective couplings -B;i,2 are given by 

B - B - iB {mN^Ai2{m%J + m^y, A21 (m^j) 

- + 2im%^ A22{m%J 

niN^ (^m^^ A2i{mlfJ + niN, Ai2{m%,j) 

B12 — B12 - 2 ±3/1 2 2 — [ — — 2 — — ? — 2"^; — • \'^-°) 

In all our results, we neglect the 1-loop corrections to the vertices W^IlNc, Zu^Na etc, 
whose absorptive parts are numerically insignificant in leptogenesis, but essential otherwise 
to ensure gauge invariance and unitarity within the PT framework [25]. 

In terms of the resummed effective couplings Bia and B'i^ and the absorptive transition 
amplitudes Aai3{s)^ the partial decay widths F^^ and their CP-conjugates P^^ are now given 
by 

= rnN^Aaa{jn\^]Bia) , P^^ = mjv, A„„(m^^; S^„) , (3.9) 

where the dependence of the absorptive transition amplitudes on Bia and has explicitly 
been indicated. Note that no summation over the individual charged leptons and light 
neutrinos running in the loop should be performed when calculating P^^ and Pj^^ using (13.30 



^Here we eliminate a typo that occurred in [9], where Ra^ in the numerator of the fraction needs be 
multipUed with —i. 
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and (13. 9p . Then, the leptonic asymmetries for each individual lepton flavour are readily 
found to be 



Ark _ 


Bla 


2 






to 


riv„ E ( 


Bla 


2 + 


B'la 





(3.10) 



r^. = E (rk + r^J , Ark = rk - rk • (s.n) 

Notice that both T^v^ and ^jy. do in general depend on the temperature T, through the 
T-dependent masses, during a second-order electroweak phase transition. More details on 
this issue will be presented in the next section. 



4 Electroweak Resonant Leptogenesis 

In this section we present the relevant Boltzmann equations (BEs) that will enable us to 
evaluate the lepton-to-photon and baryon-to-photon ratios, rj^ and r^^, during a second- 
order EWPT. In our numerical estimates, we only include the dominant collision terms 
related to the 1^2 decays and inverse decays of the heavy Majorana neutrinos A^^. We 
also neglect chemical potential contributions from the right-handed charged leptons and 
quarks [3]. A complete account of the aforementioned subdominant effects may be given 
elsewhere. 

To start with, we flrst write down the BEs that govern the photon normalised number 
densities rjj^^ and ?7al; for the heavy Majorana neutrinos No. and the left-handed leptons 
i^iL, respectively: 

(4.1) 

2 







dz 


H{Tc) 


driAL, 




dz 


H{Tc) 




2 ^k ^AL, 



(4.2) 



Although our conventions and notations follow those of [12], there are several key differences 
pertinent to our EWPT scenario that need to be stressed here. Speciflcally, we express the 
T-dependence of the BEs (14. ip and (14. 2 p in terms of the dimensionless parameter z: 

^ = ^ , (4.3) 

where Tc is the critical temperature of the EWPT to be determined below [cf. (14.60 
The parameter H{Tc) ~ 17 x T^/Mp is the Hubble constant at T = Tc, where Mp = 
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1.2 X 10^^ GeV is the Planck mass. The parameter Sjy^ denotes the branching fraction of 
the decays of the heavy Majorana neutrino into a particular lepton flavour /, i.e. — 



(r^^ + r^^)/r7VQ. Moreover, 77^^ is the equilibrium number density of the heavy neutrino 
Na, normalised to the number density of photons = 2T^/7r^: 

(4.4) 

where Kn{x) is the nth-order modified Bessel function [30]. Finally, D^^ is the T-dependent 
collision term related to the decay and inverse decay of the heavy Majorana neutrino Na- 



(4.5) 

where E^^iT) = [ IpnJ"^ + n^'N^iT) ^^"^ and = 2 is the number of hehcities of N^- 

Our next step is to include the effect of the {B + L)-violating sphalerons [4] on the 
lepton-number densities produced by the decays of Na during the EWPT. In particular, 
our interest is to implement the temperature dependence of the rate of B + L violation just 
below the critical temperature Tc, where Tc is given by [31] 

In the above, g' is the U(l)y gauge coupling and hf is the top-quark Yukawa coupling. 
We should notice that $-S mixing effects have been omitted in (14. 6p . which is a good 
approximation for scenarios with 5/A$ -C 1 as the ones to be considered here. 

A reliable estimate [32,33] of the rate of {B + L)-violating sphaleron transitions can 
be obtained for temperatures satisfying the double inequality 

Mw{T) < T < , (4.7) 

where = 5'^/4vr is the SU(2)i fine structure constant, Mw(T) = gwv{T)/2 is the 
T-dependent VT-boson mass and 

/ \ 1/2 

v{T) = v^l - —j (4.8) 

is the T-dependent VEV of the Higgs field. In detail, the rate oi B + L violation per unit 
volume is [32] 



lAiB+L) = ^Kr{UVU (^) ttg-e e-^-/^ « . (4.9) 
2 TT V 4 vr ' 
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A 


0.556 


1.612xMvi/ 


11.2 


7.6 


0.135 - 1.65 


1 



Table 1: Values of the parameters occurring in ^4-^ for ^<i>/ gl, = 0.556, which corresponds 
to a SM Higgs-boson mass of 120 GeV when 5 = 0. 

Given the double inequality (14.71) . this last expression is valid for temperatures T < T^.. 
Following the notation of [32], the parameters A/tr and A^ot that occur in (14.91) are 
functions of \<5>/gl;, VJ-ot = Svr^ and = a^T /[2 Mw{T)]. The quantity Esp is the T- 
dependent energy of the sphaleron and is determined by 

E.^^A'J^Em, (4.10) 

where A is a function of X<s>/g'^ and is 0{1), for values of phenomenological interest. The 
dependence of the parameter k on A$/(7^ has been calculated in [32,33], and the results of 
those studies are summarised in Tabled! for \q>/glj = 0.556. This value corresponds to a 
SM Higgs-boson mass Mh of 120 GeV in the vanishing limit of a ^-S mixing. 

Since the SM Higgs-boson mass is Mh 115 GeV, it can be shown [23] that the EWPT 
in the SM is not first order, but continuous from f (Tc) = to f , without bubble nucleation 
and the formation of large spatial inhomogeneities in particle densities. Therefore, we 
use the formalism developed in [12], where the {B + L)-violating sphaleron dynamics is 
described in terms of spatially independent B- and L-number densities t]b and 7]^. . More 
explicitly, the BEs of interest to us are [12]: 



drjB _ _ z Ta{b+l) 
dz ~ H{T,) 



28 f2(T) / 75 16 



(4.11) 



dVu ^ driAL, + 1 ^ (4 12) 

dz dz 3 dz ^ 

where rj^ = J2i=e,tM,TVLi is the total lepton asymmetry and 

1683 

Tacb+l) - 132T3 + 51Tv^{T) ^^^^^^^ ' ^^'^^^ 

We observe that in the limit T^(^b+l)/H{Tc) oo and for T > Tc, the conversion of 
the lepton-to-photon ratio r]L to the baryon-to-photon ratio ?7b is given by the known 
relation [34,35]: 

28 

Vb = -v^Vl- (4.14) 
51 
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Likewise, when 1 < z < 1.7 and k = 1, it is ^a{b+l)/H{Tc) ^ 1 and the baryon-to-photon 
ratio riB is then related to the total lepton-to-photon ratio rji by 

/28 16 vHT)\ r 75 t.2(r)V' „ 

For z ^ 1.7, sphaleron effects get sharply out of equilibrium and rjB freezes out. To account 
for the T-dependent {B + L)-violating sphaleron effects, our numerical estimates will be 
based on the BEs (gl]), (g^D, fHlT]) and Km . 

In the singlet Majoron model, the restoration of the global symmetry U(l); will occur 
for temperatures above a critical temperature Tj. that could in general differ from Tc of the 
SM gauge group given in (14. 6p . For example, in the absence of a doublet-singlet mixing, 
the critical temperature related to the SSB of U(l); is [27] 

r- = - 5fi . (4,16) 

A 



Consequently, the T-dependence of w{T) for T < Tj: will be analogous to v{T) in (14.81) . i.e. 



s 



^2 



1/2 



However, if m| vanishes, the singlet VEV w{T) and the doublet VEV v{T) will be related 
by an expression very analogous to (I2.10p . namely 

w{T) ^ t/v{T) . (4.18) 

As was mentioned after (I2.10p . the above relation becomes exact in a high-T expansion of 
the thermally corrected effective potential. Such an expansion is a very good approximation 
to the level of a few % for perturbatively small quartic couplings [27]. As a consequence, the 
SM gauge group and the global lepton symmetry U(l); will both break down spontaneously 
via the same second-order electroweak phase transition, with Tj: = T^. Even though the 
focus of the paper will be on this class of scenarios, we will comment on possible differences 
for models with Tj: 7^ Tc. 

If Tj = Tc, the heavy neutrino masses mTv^, the gauge-boson masses Mw,z and the 
Higgs masses Mjj s all scale with the same T-dependent factor, (1 — T'^/T^) for temper- 
atures T < Tc of our interest. Hence, the T-dependence drops out exactly in the expres- 
sion (13. 3p of the absorptive transition amplitudes Aa/3{mj^^), and likewise in the leptonic 
asymmetries and the branching fractions -BtVq- However, as can be seen from (14. 5p . the 
collision terms D^^ exhibit a non-trivial T-dependence that needs be carefully implemented 
in the BEs. 
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For our numerical estimates of the BAU, we consider the 3-generation flavour scenario 
of the RL model discussed in [11,12]. Specifically, the Majorana sector is assumed to be 
approximately SO (3) symmetric. 



rriM 



niN I3 + AM5 



(4.19) 



where AMs are small S0(3)-breaking terms that are of order rn'^m£)/m]\f as these are 
naturally expected from (12.201) . Plugging (14.191) into (I2.20p . we find that, to leading order 
in AM5, the heavy neutrino mass matrix deviates from mjv I3 by an amount 



= AMs + i'm^D^D + "^d"^d) • 

2mN ^ ' 



(4.20) 



It is interesting to observe that possible renormalisation-group (RG) running effects from 
a high-energy scale Mx, e.g. GUT scale, down to will induce a negative contribution 
to [18], i.e. 



(5m^)RG 



(^m\jmD + ^E"^!) ) 111 



M 



X 



rriN , 



(4.21) 



For Mx = Mgut ~ 10 and wln = 80-150 GeV, the RG-induced terms are typically 
smaller by a factor ~ 0.1-0.4 with respect to the tree-level contribution given in (I4.20p . 
Thus, the inclusion of the RG effects are not going to affect the results of our analysis in 
a substantial manner. 

As was mentioned already, the SO (3) symmetry is broken by the Dirac mass terms 
(^z))io, which in our case possess an approximate U(l)-symmetric flavour pattern [11]: 



rriD 



V 

V2 



( ae-^^l^ 
\ ce^*''/^ ce^''/'^ 



(4.22) 



where the 3-by-3 matrix ^m^i. 



5mn = —= 



V 

V2 



\ 




\ Er / 



(4.23) 



violates the U(l) symmetry by small terms of order of the electron mass rrie. Instead, the 
U(l)-symmetric Yukawa couplings a and b can be as large as the r-lepton Yukawa coupling 
rrir/v, i.e. of order 10^^-10^^. For successful RL, it was found [11,12] that the parameter 
c needs to be taken of the order of the electron Yukawa coupling rrie/v. It is important to 
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Higgs 
Sector 



tan/3 



Mh [GeV] 



Ms [GeV] 



0.238 



30 



1/V2 



121 



29 



Neutrino 
Sector 



(^m^)ii 
rriN 



((^m^)i2 



(<5m^)i3 

TUN 



TTT-AT 



((^m^)23 

mAT 



(^m^)33 



10" 



-10- 



-4 X 10 



-10 



4 X 10" 



(6.8 - 0.6i) 
xlO-9 



5.2 X 10" 



3 
5DU 



57 



2 X 10 



-7 



1563 
250000 



39 
30000 



100 GcV 



100 GeV 



100 GcV 



100 GcV 



147 
128000 

/ Amjv ' 

/ 100 GeV 



Table 2: Complete set of the theoretical parameters used for the singlet Majoron model, 
where ArriN = 2((5m^)23 + z[(5m^)33 - (5111^)22] ■ 

stress here that the approximate flavour symmetries S0(3) and U(l) ensure the stabihty 
of the light- and heavy-neutrino sector under loop corrections [11,13,36]. 

For our numerical analysis, we fully specify in Table [2] the values of the theoretical 
parameters for the Higgs and neutrino sectors. The only parameter that we allow to vary 
is the heavy Majorana mass scale m^. For 50 GeV < m^v ~ 200 GeV, the choice of 
parameters in Table [2] leads to an inverted hierarchical light-neutrino spectrum with the 
following squared mass differences and mixing angles: 

ml^ - m\ = (7.5-7.7) x 10"^ eV^ , ml^ - ml^ = 2.44 x 10"^ eV^ , 

sin2^i2 = 0.362, sin2^23 = 0.341, sin^ ^13 = 0.047 (4.24) 

and mj,3 = 0. The spectrum is compatible with the light-neutrino data at the 3a confidence 
level (CL) [37]. 

In Fig. [2] we present numerical estimates of the lepton-fiavour asymmetries rji^^^ 
and the baryon asymmetry rjs &s functions oi z = T^/T, for a typical electroweak RL 
scenario with niN = 100 GeV. As initial conditions at T = Tc ~ 133 GeV, we take 
?7™^ = 1 for the heavy neutrino number densities and vanishing lepton-to-photon and 
baryon-to-photon ratios, i.e. rjf^^ ^ ^ = ^^^^1 = 0. The thermal in-equilibrium condition 
77™^ = 1 is expected, since the heavy neutrinos A'^1,2,3 have no chiral masses when T > 
Tc and get rapidly thermalised by the sizeable light-to-heavy neutrino Yukawa couplings 
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Figure 2: Numerical estimates of tjb (solid), rj^^ (dash- dotted), rji^ = rj^^ (dotted) and 
TjL (dashed) as functions of z = T^jT , for a model with niN = 100 GeV, and rjf}^ = 1. 
The model parameters are given in Table The horizontal grey line corresponds to the 
observed baryon-to-photon ratio rj'^'^ = 1.65 x 10^^, after evolving the latter back to the 
higher temperature T = Tc/10. 

can be seen from Fig. [2], a net baryon asymmetry r^^ is generated 
by a non-zero r-lepton asymmetry rj^^ . Tliis LT--excess is created before sphalerons sharply 
freeze out, i.e. for temperatures T > Tgph ~ 78 GeV {z < 1.7). Consequently, in the thermal 
evolution of the Universe, there is a sufficiently long interval 78 GeV ^ T < 133 GeV, 
where a leptonic asymmetry can be converted into the observed BAU for our scenarios 
with spontaneous lepton-number violation at the electroweak scale. 

Figure [3] exhibits the dependence of the baryon-to-photon ratio tjb on z = T^/T 
for different values of the heavy Majorana mass scale rriN. We notice that the lighter 
the heavy neutrinos are, the smaller the created baryon asymmetry is. For example, for 
heavy-neutrino masses m^v ~ 80 (50) GeV, tjb falls short almost by one order (two orders) 
of magnitude with respect to the observed BAU rj]^^. This is a generic feature of our 
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Figure 3: Numerical estimates of rjB versus z = T^/T for = 120 GeV (dashed), 
100 GeV (solid), 80 GeV (dash- dotted), 50 GeV (dotted). The meaning of the horizontal 
grey line is the same as in Fig. [H 

electroweak RL scenarios based on large wash-out effects due to the relatively large Dirac- 
neutrino Yukawa couplings {mD)ia/v. If the heavy neutrinos have masses rriM < 90 GeV, 
their number densities will start decreasing for T < niN, potentially creating a net lepton 
asymmetry that can be converted into rj^^. However, this should happen above the freeze- 
out temperature Tgph ~ 78 GeV of sphalerons. Thus, successful electroweak RL requires 
that rriN > T^pJ^ 

■'■Recently, a different leptogenesis scenario ■with to^v <C TLph was studied in [38], ■where the BAU is 
generated by sterile-neutrino oscillations. Such a realisation relies on the assumption that the oscillating 
sterile neutrinos start evolving from a coherent state and retain their coherent nature within the thermal 
plasma of the expanding Universe. In the singlet Majoron model ■we have been studying here ho^wever, 
t-channel 2^2 scattering processes, such as JJ ^ ^aR'^aR, that occur before the EWPT (T > Tc) are 
strong, ■with rates 0[p^g^T / (Stt)] 3> H{T), for Higgs-singlet Yuka^wa couplings paa ^ 1- They can therefore 
lead to rapid thermalization and loss of coherence of the massless right-handed neutrinos. Shortly after 
the EWPT, for z = TJT > 1.1, it is Tn^,.,IH ~ lO^-lO^o and TnJH 1-10, which again gives rise to an 
almost instant thermalization of all the heavy neutrino mass eigenstates A^i,2,3- 
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Finally, it is important to comment on the last condition > Tgph for scenarios 
with 7^ Tc. This condition will still be valid, as long as Tj: > Tgph- However, for scenarios 
with Tj: < Tc, the predicted BAU rjB will sensitively depend on the initial values vTef^r 
rf^aXT = Tc- Instead, if T^' ^ Tc and > 90 GeV, the predictions for the BAU will 
remain almost unaffected, even if r/^ ~ 7]°^^ at T = 10 Tc [12]. 

5 Astrophysical and Phenomenological Implications 

It is interesting to discuss the implications of the singlet Majoron model for astrophysics 
and low-energy phenomenology. To quantify the effects of heavy Majorana neutrinos, we 
define the new-physics parameters 

^iv = — Bif^Biik = Bi^Bi>a , (5.1) 

where I, I' = e, fi, r. Evidently, in the absence of light-to-heavy neutrino mixings, the 
parameters Qur vanish. LEP and low-energy electroweak data put severe limits on the 
diagonal parameters Qu [39]: 

fice < 0.012, n^^ < 0.0096, n^r < o.oie, (5.2) 

at the 90% CL. On the other hand, lepton-fiavour-violating (LEV) decays, such as /i ^ 
67 [40], fi eee, t ^ e'j, t ^ eee, /z — >• e conversion in nuclei [41,42] and Z 
IV [43], constrain the off-diagonal parameters Vtui^ with I ^ V . The derived constraints 
strongly depend on the heavy neutrino masses ^«Ar„ and the size of the Dirac masses {mj:,)ia. 
However, for models relevant to leptogenesis, with {mD)ia ^ Mw [41], we obtain the 
following limits: 

I^cmI ^ 0.0001, l^crl ^ 0.02, \n^,r\ ^ 0.02, (5.3) 

including the recent BaBar data on LEV r decays [44]. 

The predictions for LEV decays in models of resonant leptogenesis has been exten- 
sively discussed in [12]. Since our results obtained in Section H] agree well with this earlier 
analysis, we will not repeat the details of this study here. Here, we only reiterate the fact 
that successful electroweak RL requires that 171^ > 100 GeV. This latter constraint gives 
rise to the following upper limits: 

fiee ^ 2.2 X 10"^ , I^cmI ^ 8.3 X 10"^ , n^f, < 3.1 X 10"^ , (5.4) 

whereas all remaining parameters Qw are O{10^^) and so unobservably small. All these 
limits are deduced by using the model parameters of Table [2l 
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Figure 4: Loop-induced couplings of the Majoron to charged leptons I, I' and quarks q. 

In the singlet Majoron model under study, there are additional LFV decays for the 
muon and the tau-lepton that involve the Majoron, i.e. /i — > Je, r Je and r Jfi. As 
shown in Fig. |U these LFV decays are induced by heavy Majorana neutrinos at the 1-loop 
level. Detailed analytic expressions for the loop-induced couphngs JW and Jqq, where q is 
a quark, may be found in [21]. To leading order in Qui, the prediction for the LFV decay 
/- ^ I'- J is 

m^i'j) - ^r-^\ = '-^ (i ^ (5.5) 

where \n = 'm\i/M'y^. For Aa? = 1, the prediction for the observable R{1 I' J) takes on 
the simpler form: 



3«u, ,2 in |2 



= ^4|0„f ^. (5.6) 

The requirement for successful electroweak RL, i.e. > 100 GeV, gets translated into 
the following upper bounds: 

R{fi eJ) < 2.7 X 10"^ , R{t eJ) < 4.6 x 10"^^ , R{t ^J) < 6.7 x 10^^^ . 

(5.7) 

On the experimental side, however, the following upper limits are quoted: 



It is interesting to remark that the predicted value for R{fi — > eJ) is close to the present 
experimental sensitivity, whereas the other decay modes turn out to be very suppressed for 
the given RL model with inverted light-neutrino hierarchy. Had we chosen a model with 
normal hierarchy, the decay rates R{t — > eJ) and R{t — > fj,J) would have been enhanced 



R{fi - 


eJ) 


< 


2.6 X lo~^ 


at 90% CL [45] 


R{t - 


eJ) 


< 


1.5 X lo~^ 


at 95% CL [46] 


R{t - 




< 


2.6 X 10"% 


at 95% CL [46] 
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by a factor ~ 10^, but they will still be rather small (9(10 ^) to be observed; the predictions 
generally lie 4 orders of magnitude below the current experimental upper bounds. 

Useful constraints on the parameters of the theory are obtained from astrophysics as 
well [47]. Specifically, observational evidence of cooling rates of white dwarfs implies that 
the interaction of the Majoron to electrons, gjeeJei'y^e, should be sufficiently weak and the 
coupling gjee must obey the approximate upper bound [48]: 

\9jee\ < 10-^' . (5.9) 

The above limit gets further consolidated by considerations of the helium ignition process 
in red giants, leading to the excluded range: 3 x 10"^^ ^ \gjee\ ~ 6 x 10~^. To leading 
order in flu, the loop-induced coupling gjee, is given by [21]: 



djee = —rp. -Tr~ '^P 



1 - Aa/ \ 1 - Aa/ / 2 , ^ 



\N \ ^ — 

If Aat ^ 1, the expression for the coupling gjee simplifies to 



(5.10) 



whilst for Aat = 1 gj^e becomes 

Given the limits (15.41) for successful RL, we can estimate that 

(?jee ^ -3.3 X 10-1^ , (5.13) 

which passes comfortably the astrophysical constraint given in (15. 9p . 

Useful astrophysical constraints may also be obtained from considerations of the cool- 
ing rate of neutron stars [47]. Neutron stars will loose energy by Majoron emission through 
the interaction: g,jj\fj\f J Wi'-f^N', where A/" is a nucleon, specifically a neutron. The obser- 
vational limit on gjj\fj\f is [49] 

g^M ^ 10-' . (5.14) 
On the other hand, the theoretical prediction for gjqq at the quark level is 

gjqq = ^^2^" "^^ (^^'^ ~^ ~^ ^tt) ■ (5.15) 

From naive dimensional analysis arguments, one expects that gjj\fj\f ~ {mj^/mg)gjgg. In 
this way, one may estimate that 

gjNM ~ 7 X 10-1° , (5.16) 



21 



after taking into consideration the limits stated in (15.41) . 

Cosmic microwave background (CMB) data and BBN put stringent limits on the max- 
imum number of weakly-interacting relativistic degrees of freedom, such as light neutrinos 
and Majorons [50,51]. In particular, the allowed range obtained for the effective number 
Ny of left-handed neutrino species is = 2.70l5';32 at the 68% CL [51]. The upper bound 
on may naively be translated into an upper limit on AN^ = N^, — 3 = 0.61 of extra 
effective neutrino species beyond the 3 SM left-handed neutrinos. The singlet Majoron 
contributes AN,^ = (| x j)^^^ ~ 0.474, if its freeze-out or decoupling temperature Tj is 
equal to the corresponding one T^, of the neutrinos. Although this result does not pose 
by itself a serious limitation on the singlet Majoron model, it can be estimated, however, 
that Tj ^ Tj, ^ 1 MeV and the contribution of J to AN^ becomes even more suppressed. 
Specifically, the freeze-out temperature Tj is determined when the annihilation rate of 
Majorons through the process J J vv becomes smaller than the Hubble expansion rate 
H{T) of the Universe. The annihilation process J J vv is mediated by the H and S 
bosons in the s-channel and by the heavy neutrinos A^^i 2,3 in the t-channel. Considering 
the latter reactions only, one may naively estimate that 



A similar value for Tj/T^ is obtained if the S*, if -boson exchange processes are used for 
the model parameters of Table [2l Thus, the freeze-out temperature Tj lies in the range 
0.1-1 GeV, namely about the quark-hadron deconfinement phase. In this epoch of the 
Universe, the effective number of relativistic degrees of freedom is 5'*(Tj) ^ 66. Then, the 
actual contribution of the Majoron to AA^",^ is reduced with respect to the Tj = T^, case 
by a factor {g^{T^) / g^{Tj)f'^ ^ 0.016 to the value AA^^ ^ 0.008, which is far below the 
present and future observational sensitivity [51]. 

Finally, singlet Majorons J and singlet scalars S may also give rise to interesting 
collider phenomenology [52] through the singlet-doublet mixing parameter 5 in the scalar 
potential (12.11) . However, since 5 -C A$ (cf. Table [2]), the singlet Majoron scenario under 
study predicts a rather small mixing angle sq ~ —0.1. The production cross section of S, 
via the process e+e~ — > ZS, is then suppressed with respect to the SM one by a factor 
~ 0.01. Moreover, the so-produced Higgs singlets may decay quasi-invisibly into a pair 
of Majorons J, which makes difficult to fully rule out such a scenario by LEP2 data or at 
the LHC. Future high-energy e+e~ colliders of higher luminosity will severely constrain the 
allowed parameter space of this singlet Majoron model. 



7j 




(5.17) 
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6 Conclusions 



The origin of CP violation in nature still remains an open physics question. If CP violation 
originates from the SSB of the SM gauge group, the original scenario [1] of GUT-scale 
leptogenesis will be excluded. Similar will be the fate of all high-scale leptogenesis models, 
if the source of lepton-number violation is due to the SSB of a global U(l); symmetry at the 
electroweak scale. In this paper we have shown how resonant leptogenesis at the EWPT 
constitutes a realistic alternative for successful baryogenesis in models with spontaneous 
lepton-number violation. Specifically, we have considered a minimal extension of the SM, 
the singlet Majoron model, which includes right-handed neutrinos and a complex singlet 
field that carries a non-zero lepton number. Depending on the form of the scalar potential, 
the lepton number can get broken spontaneously through the VEV of the SM Higgs doublet. 
Taking into consideration the Boltzmann dynamics of sphaleron effects, we have analysed 
the BAU for different values of the Majorana mass scale m^r within the context of a 
benchmark scenario whose model parameters are given in Table El The generic constraint 
from having successful electroweak RL is that rriN ^ Tsph, where Tgph ~ 78 GeV is the 
freeze-out temperature of the sphalerons. 

The singlet Majoron model predicts a massless Goldstone particle, the Majoron J. 
The Majoron can be produced via the LFV decays, /i Je, t ^ Jji and r — > Je. 
Considering the constraints from successful electroweak RL and the astrophysical limits 
derived from the cooling rate of neutron stars, we have found that the decay mode ^ ^ Je 
is the most promising channel, with sizeable branching fraction that can be looked for in 
the next-round low-energy experiments. 

The predictions obtained for the BAU in this study are limited by the approximations 
that are inherent in the calculation of the non-perturbative sphaleron dynamics. The 
predicted values should be regarded as order-of-magnitude estimates, since the {B + L)- 
violating sphaleron transitions crucially depend on the parameter n that varies by a factor 
of 10 or so. It would therefore be very valuable to go beyond the current approximation 
methods and improve the computation of the out-of-equilibrium sphaleron dynamics during 
a second-order electroweak phase transition. 
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